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ABSTRACT 

The  J-Integral  is  an  elastic -plastic  fracture  criterion,  which  permits  measurement  of  the 
fracture  toughness  of  a  specimen  that  has  been  fractured  after  general  yielding.  However,  an 
understanding  of  the  ratio  of  plastic  zone  size  (radius)  to  the  crack  tip  blunting  (stretch  zone) 
required  to  determine  the  upper  limit  of  temperature  relative  to  full  size  transition  curves  where 
elastic  plastic  fracture  becomes  invalid.  This  study  endeavours  to  acquire  this  ratio  using  finite 
element  techniques. 

The  development  of  the  plastic  zone  in  dynamic  tear  (DT)  specimens  and  a  non-  standard 
three  point  bending  fracture  test  specimen  used  to  measure  fracture  properties  was  the  main 
focus  of  the  study.  The  ABAQUS  finite  element  software  was  used  to  model  the  elastic-plastic 
behaviour  of  the  specimens.  For  the  DT  specimen,  a  crack  was  induced  by  pressing  the  notch, 
followed  by  fatigue  cracking  at  a  limit  load  level  of  40%  of  the  specimen  limit  load,  whereas,  the 
crack  shape  for  the  non-standard  specimen  was  a  fatigue  crack  defined  at  approximately  30%  of 
the  limit  load.  The  shapes  of  these  cracks  were  adequately  modelled  in  the  finite  element 
analysis.  The  specimens  were  made  of  350WT  steel  and  304  stainless  steel  materials.  The 
specimens  were  loaded  until  fixed  amounts  of  permanent  deformation  were  recorded.  Results 
were  obtained  in  the  form  of  plots,  showing  the  progression  of  the  plastic  zone  around  the  crack 
tip.  For  each  case,  the  results  provide  the  following:  mid  point  plastic  deflection,  stretch  zone 
width  and  plastic  zone  radius.  The  finite  element  results  obtained  were  compared  to 
experimental  elastic  plastic  testing  where  available,  and  reasonably  accurate  agreement  was 
achieved. 

RESUME 

L’integrale-J  est  un  critere  de  rupture  Elastique-plastique  qui  permet  de  mesurer  la 
tenacity  d’un  specimen  qui  a  subi  une  rupture  suite  a  une  deformation  generale.  D  est  toutefois 
necessaire  de  comprendre  le  rapport  de  la  dimension  de  la  phase  plastique  (le  rayon)  au  degre 
d’emoussement  de  la  tete  de  fissure  (la  zone  d’allongement)  afin  d’Etablir  la  limite  superieure  de 
temperature,  par  rapport  aux  courbes  de  transition  a  l’Echelle,  a  laquelle  la  rupture  elastique- 
plastique  n’est  plus  valable.  La  presente  etude  s’efforce  de  determiner  ce  rapport  a  l’aide  de 
1’ analyse  par  elements  finis. 

L’ etude  mettait  1’ accent  sur  la  croissance  de  la  phase  plastique  dans  des  Eprouvettes  ay  ant 
subi  l’essai  dynamique  de  dechirement  (eprouvettes  DT),  ainsi  que  dans  une  eprouvette  non 
standard  ayant  subi  un  essai  de  rupture  avec  pliage  en  trois  points  et  utilisee  pour  mesurer  les 
proprietes  de  fracture.  La  modelisation  du  comportement  Elastique-plastique  des  eprouvettes  a 
ete  effectuee  a  l’aide  du  logiciel  d’analyse  par  elements  finis  ABAQUS.  Dans  le  cas  de 
1’eprouvette  DT,  une  fissure  a  ete  causee  par  la  pression  exercEe  sur  l’entaille,  suivie  d’une 
fissuration  de  fatigue  a  un  niveau  Equivalent  &  40  %  de  la  charge  limite  de  1’ eprouvette,  alors  que 
la  forme  de  la  fissure  dans  l’eprouvette  non  standard  correspondait  a  une  fissure  de  fatigue 
equivalant  a  environ  30  %  de  la  charge  limite.  L’ analyse  par  elements  finis  a  permis  de  bien 
modeliser  la  forme  de  ces  fissures.  Les  eprouvettes  etaient  composees  d’acier  350WT  et  d’acier 
inoxydable  304.  Elies  etaient  mises  a  1’ essai  jusqu’a  l’obtention  de  valeurs  constantes  de 
deformation  permanente.  Les  resultats  obtenus  montrent,  sous  forme  graphique,  la  propagation 
de  la  phase  plastique  dans  la  tete  de  fissure.  Les  resultats  foumissent,  pour  chaque  cas,  le 
flechissement  plastique  au  point  milieu,  la  largeur  de  la  zone  d’allongement  et  le  rayon  de  la 
phase  plastique.  Les  resultats  de  1’ analyse  par  elements  finis  ont  ete  compares  aux  donnees 
experimentales  disponibles  pour  des  essais  de  deformation  elastique-plastique  et  la  concordance 
Etait  raisonnablement  precise. 
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1.1 


1.  INTRODUCTION 

1.1  Background 

The  J-integral  is  one  of  several  criteria  used  to  determine  the  fracture  toughness  of 
materials.  This  is  a  criterion  that  is  suitable  for  measurement  of  fracture  in  the  elastic-plastic 
regime  of  the  full  transition  curve  (Figure  1.1).  It  permits  measurements  of  fracture  toughness  of 
specimens  that  have  been  fractured  after  general  yielding,  and  assumes  that  onset  of  fracture  is  due 
to  a  critical  level  of  stress  or  strain  being  achieved  at  or  near  the  crack  tip  [1,2].  The  J-integral  is 
measured  using  three-point  bend  specimens,  such  as  the  dynamic  tear  (DT)  test  specimens.  An 
understanding  of  the  ratio  of  the  plastic  zone  size  (radius)  to  the  crack  tip  blunting  (stretch  zone 
width)  is  required  to  determine  the  upper  limit  of  temperature  relative  to  the  full  size  transition 
curves  where  elastic-plastic  fracture  (J-integral)  becomes  invalid. 

Figure  1 .2  shows  a  typical  DT  specimen  and  Figure  1 .3  shows  idealizations  of  the  plastic 
zone  around  the  crack  tip.  As  shown  in  Figure  1 .3,  the  plastic  zone  radius  is  the  maximum  extent 
of  the  plastic  zone,  and  the  stretch  zone  width  (SZW)  is  the  distance  from  the  crack  tip  to  the 
interface  of  the  plastic  zone  with  the  side  of  the  crack.  The  shear  lip  size  is  given  by  the  distance 
from  the  yield  surface  to  the  original  surface  on  the  side  of  the  specimen.  The  standard  validity 
criterion  for  stipulates  that 


a,b,B  =  50J/<yJlow  (1.1) 

This  suggests  that  B  be  greater  than  50J/aflow  But  blunting  is  approximately  J/(2aflow).  Therefore 
the  thickness  B  should  be  greater  than  100  times  the  critical  blunting  for  a  valid  test.  Thus,  if  the 
ratio  of  plastic  zone  size  to  blunting  size  is  known,  then  the  following  relations  will  hold  for  a 
valid  fracture  test: 


B>\  00SZW 

B  >1  OOr^,  / P  (1.2) 

B  >  100s7(3  (or  s  <  5(3/100) 
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where  ry  =  plastic  radius,  P  =  ry/SZW,  S  is  the  shear  tip  size  and  the  flow  stress  oflow  is  given  by 
°flow  =  (oy+ou)/2,  where  ov  is  the  yield  stress  and  ou  the  ultimate  stress.  Accordingly,  in  the  present 
study,  the  key  is  to  identify  the  ratio  of  plastic  zone  size  to  stretch  zone  width.  Four  methods  of 
verification  are  available: 

(i)  Analytical  (after  Rice[3],  who  suggested  the  ratio  was  as  low  as  5  and  as  high  as  600; 

(ii)  Visual  (after  experiments  on  stainless  steel  at  DREA); 

(iii)  Experience  in  valid  J  testing  (after  DREA  investigations);  and 

(iv)  By  finite  element  (looking  at  plastic  deformation  of  the  center  of  the  bar,  stretch  zone  or 
blunting,  plastic  zone,  and  the  parameter  s  which  actually  includes  the  shear  lip  and  lateral 
contraction). 

1.2  Objectives  and  Scope 

The  focus  of  this  study  is  to  utilize  the  finite  element  methodology  to  investigate  the 
plastic  zone  development  in  three  point  bend  fracture  specimens.  The  aim  of  the  study  is  to 
provide  an  understanding  of  the  ratio  of  the  plastic  zone  size  (radius)  to  crack  tip  blunting  (stretch 
zone  width),  in  order  to  be  able  to  predict  the  upper  limit  where  elastic  plastic  fracture  becomes 
invalid.  Of  interest  is  the  investigation  of  the  plastic  zone  development  in  standard  dynamic  tear 
(DT)  test  specimens  with  350WT  or  304  stainless  steel  materials;  and  a  non-standard  (NS)  three- 
point  bend  specimen  made  of  304  stainless  steel  material. 

Chapter  2  provides  detailed  descriptions  of  the  specimen  configurations  considered  in 
the  study,  and  a  description  of  the  crack  tip  shape.  The  finite  element  methodology  utilized  in  the 
investigation  is  presented  in  Chapter  3,  which  provides  details  of  the  finite  element  meshes, 
materials  constitutive  models,  boundary  and  loading  conditions,  and  the  solution  methodology. 
The  finite  element  results  are  presented  in  Chapter  4,  with  discussions  of  the  results.  Finally, 
Chapter  5  provides  a  summary  of  the  study  and  the  conclusions  reached. 
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Figure  1.1:  Typical  Transition  Curve 
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Figure  1.2:  Configuration  of  DT  Specimen 
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2.1 


2.  PROBLEM  DESCRIPTION 

2.1  Specimen  Configurations 

A  schematic  representation  of  the  test  specimen  was  presented  in  Figure  1 .2.  In  the 
present  study,  three  test  specimens  are  analyzed.  These  include  the  following: 

(i)  DT  specimen  with  350WT  steel,  designated  as  DT-350WT; 

(ii)  DT  specimen  with  304  stainless  steel,  designated  as  DT-304SS;  and 

(iii)  Non-standard  specimen  with  304  stainless  steel,  designated  as  NS-304SS. 

The  DT  specimens  (DT-350WT  and  DT-304SS)  had  the  following  geometric  properties: 

L  =  181  mm; 

S  =  164  mm; 

a  =  12  mm  plus  .010  inches  pressing  followed  by  fatigue  extension  of  1 3  to  20  mm; 

W  =  ~41  mm;  and 

B  =  8  mm. 

The  non-standard  specimen  had  the  following  configuration: 

L  =  125  mm; 

S  =  100  mm; 

a  =  12.5  mm  fatigue  crack  produced  by  rigorous  process  (to  be  included  in  FE  model); 

W  =  25  mm;  and 

B  =  12.5  mm 

Where 

L  =  length  of  beam; 

S  =  span  of  beam  between  supports; 

W  =  width  of  beam; 

B  =  thickness  of  beam;  and 

a  =  notch  size  (fatigue  sharpened) 

2.2  Crack  Tip  Shape 

The  crack  tip  shape  is  that  of  a  fatigue  crack  initially  and  a  blunted  fatigue  crack 
thereafter. 
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For  the  DT  specimen  the  notch  is  first  machined  with  an  included  angle  of  60°  followed 
by  a  pressing  process  that  extends  the  notch  0.254  mm  (0.010)  with  an  included  angle  of  40°  and 
a  resulting  crack  tip  radius  of  0.254  mm  (0.001).  This  is  followed  by  fatigue  cracking  at  a  load 
level  of  no  more  than  40%  of  the  current  specimen  limit  load.  This  produces  a  sharp  fatigue  crack 
of  definable  dimensions,  which  we  would  like  to  use  in  the  FE  simulation.  This  pre-cracking 
extends  the  pressed  notch  by  at  least  1  mm. 

The  crack  tip  shape  for  the  non-standard  specimen  is  that  of  a  well  established  fatigue 
crack  whose  maximum  loading  has  been  defined  by  careful  laboratory  fatiguing  (maximum  load 
approximately  30%  of  the  limit  load).  The  shape  of  the  resulting  crack  tip  and  the  configuration 
of  the  crack  along  its  entire  length  were  acquired  by  sectioning.  Changes  in  shape  of  the  crack  tip 
as  the  FE  model  is  run  from  one  plastic  boundary  condition  to  another  are  compared  to  detailed 
laboratory  metallo graphic  measurements.  Detailed  descriptions  of  the  precracking  process  is 
available  in  Reference  [4]. 

2.3  Requirements 

It  was  required  to  compute  the  progression  of  the  plastic  zone  around  the  crack  tip  under 
quasi-static  loading  conditions,  for  the  three  specimens.  For  each  case,  the  following  response 
parameters  were  required: 

(i)  Mid  point  plastic  deflections  relative  to  S  and  L  mentioned  above; 

(ii)  Stretch  zone  width  (i.e.  distance  from  crack  tip  to  interface  of  yield  zone  with  side  of 
crack);  and 

(iii)  Plastic  zone  radius  (i.e.  locus  and  maximum  extent  of  plastic  zone). 
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3.  FINITE  ELEMENT  MODEL 

3.1  Finite  Element  Approach 

The  prediction  of  the  response  behaviour  of  the  DT  specimen  requires  a  tool  that 
accounts  for  large  geometry  change  and  plasticity  effects  in  the  crack  tip  region.  There  are  several 
commercial  FEA  tools  that  can  be  employed  or  adapted  for  calculating  the  stresses  and  strains 
needed.  In  this  study  the  ABAQUS  finite  element  software  [5]  was  used  to  model  the  elastic- 
plastic  behaviour  of  the  DT  specimen.  The  software  has  a  wide  range  of  non-linear  materials 
models  and  can  account  for  large  strains  and  displacements.  Moreover,  it  has  been  successfully 
used  for  the  analysis  of  similar  problems  [6-8].  However,  the  HyperMesh  [9]  general-purpose  pre- 
and  post-processing  program  was  used  for  model  generation  and  results  processing.  Details  of  the 
finite  element  model  are  provided  below. 

3.2  Finite  Element  Meshes 

As  stated  above,  the  HyperMesh  code  was  used  to  generate  the  finite  element  models  of 
the  DT  and  non-standard  specimens,  which  were  then  translated  to  ABAQUS  input  files.  The 
exact  shapes  of  the  crack  were  obtained  from  DREA  to  enable  accurate  modelling  of  the  structural 
configuration.  Figure  3.1  shows  the  shape  of  the  fatigue  crack  in  the  DT  specimen.  This  shape  had 
to  be  modelled  accurately  to  provide  meaningful  results.  Detailed  2-D  finite  element  models  of 
the  specimen  were  developed.  In  order  to  reduce  the  problem  size  only  one-half  of  the  structure 
was  modelled.  The  finite  element  models  of  the  DT  and  non-standard  (NS)  specimens  are  as 
shown  in  Figures  3.2  and  3.3,  respectively.  Plane  strain  elements  were  used  to  model  the  structure. 
These  elements  allow  for  the  treatment  of  large  displacements,  finite  strains  and  plasticity,  which 
are  expected  to  occur  in  the  specimen. 

3.3  Material  Models 

An  incremental  rate  independent  plasticity  theory  available  in  the  ABAQUS  finite 
element  program  [5]  was  used  for  the  material  constitutive  model.  This  standard  model  for 
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3.2 


plasticity  is  summarized  here  for  completeness.  The  total  multi-axial  strain  state  eu,  expressed 
in  terms  of  elastic  and  plastic  components,  was 


:*  4- 

if  +  SU 


(3.1) 


The  total  logarithmic  uniaxial  strain,  e,  consistent  with  the  integration  of  the  rate  of  deformation 
tensor  for  a  multiaxial  strain  state,  is  decomposed  as 


€  =  €  +  S' 


The  yield  function  f  is  related  to  uniaxial  tension  by 


(3.2) 


(3.3) 


where  Ty  and  t  are  the  multiaxial  Kirchoff  and  uniaxial  stress  states,  respectively.  The  associated 
flow  rule,  dissipation  equivalence  condition  and  a  consistency  condition  govern  plastic  strain 
increments. 


The  specific  form  of  f  is  the  von  Mises  yield  function 

/,  \  1/2 

/(’»)- (|(v*)J  '  <3-3> 

where  Sy  is  the  deviatoric  component  of  stress. 

The  Kirchoff  stress  and  logarithmic  strain  measures  are  employed  because  of  advantages 
gained  in  computational  implementation.  The  Kirchoff  stress  tensor,  Xy,  is  approximately  equal 
to  the  more  physically  motivated  Cauchy  stress  tensor,  Oy,  for  deformations  involving  only  small 
changes  in  volume  [7,8].  This  condition  is  implicit  in  the  current  elastic-plastic  analysis.  The 
uniaxial  Cauchy  stress-logarithmic  strain  constitutive  response  of  the  material  are  formally 
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input,  in  multilinear  form,  as  Cauchy  stress  and  logarithmic  plastic  strain  paris  for  the  ABAQUS 
program.  The  constitutive  parameters  for  the  350WT  and  304SS  materials  are  described  below. 

3.3.1  350WT  Steel 

The  material  for  the  first  DT  specimen  was  350WT  steel.  The  engineering  stress-strain 
curve  for  this  material  was  obtained  from  DREA.  This  curve  terminated  at  a  strain  of  25%. 
However,  from  experimental  investigations  by  DREA  [10],  strains  as  high  as  250%  were 
observed  around  the  crack  tip  of  the  specimen.  A  method  for  defining  the  stress-strain  behaviour 
beyond  25%  strain  was  therefore  required.  Two  possible  stress-strain  curves,  as  illustrated  in 
Figure  3.4,  were  considered.  In  the  first  model  (Figure  3.4a),  the  stress  was  assumed  to  be 
constant  at  the  maximum  stress  level  from  25%  to  250%  strain.  This  case  is  referred  to  as  the 
flat  curve  model.  In  the  second  constitutive  model  (Figure  3.4  (b)),  the  stress  was  assumed  to 
rise  from  the  maximum  uniaxial  test  stress  of  about  600  MPa  to  1100  MPa  from  25%  to  250% 
strain.  The  maximum  stress  value  was  determined  from  experimental  observation  and  intuition 
[10].  This  second  case  is  referred  to  as  the  rising  curve  model. 

3.3.2  304  Stainless  Steel 

The  second  DT  specimen  and  the  non-standard  specimens  were  made  of  304  stainless 
steel  material.  The  stress  strain  curve  for  the  material  was  obtained  from  the  literature  [11],  This 
curve  terminated  at  a  strain  of  about  80%.  However,  from  experimental  investigations  by  DREA 
[5],  strains  as  high  as  1000%  were  observed  around  the  crack  tip  of  the  specimen.  In  this  case 
only  the  rising  curve  model  was  used  beyond  the  80%  strain  level.  The  possible  stress-strain 
curve  is  illustrated  in  Figure  3.5.  In  this  constitutive  model  the  stress  was  assumed  to  rise  from 
the  maximum  uni-axial  test  stress  of  about  600  MPa  to  1000  MPa  from  80%  to  1000%  strain. 
The  maximum  stress  value  was  determined  from  experimental  observation  and  intuition  [10], 
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Boundary  Conditions  and  Loading 


The  following  boundary  conditions  were  applied: 

At  the  support:  v  =  0 

Along  the  center  line:  u  =  0; 

Where,  u,v  are  the  displacements  components  in  the  longitudinal  and  transverse  directions.  The 
load  was  applied  as  a  concentrated  load  at  the  top  middle  point  of  the  beam. 
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Figure  3.2:  FE  Mesh  of  DT  Specimens 
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4.  RESULTS  AND  DISCUSSIONS 

This  chapter  discusses  the  results  of  the  study.  The  finite  element  results  obtained  for 
the  three  specimen  configurations  are  discussed  in  Sections  4.1  to  4.3,  respectively.  For  each 
specimen,  results  of  the  midpoint  displacement,  plastic  zone  development  and  stretch  zone  width 
are  provided. 

4.1  DT  Specimen  with  350  WT  Steel  Material 

Midpoint  Displacements 

The  predicted  load-displacement  response  of  the  midpoint  of  the  beam  is  shown  in 
Figure  4.1.  In  the  figure,  DL  refers  to  the  displacement  of  the  top  midpoint  with  respect  to  the 
top  left  end  point,  and  DS  refers  to  the  bottom  midpoint  with  respect  to  the  support  point.  The 
results  obtained  from  using  the  flat  and  rising  stress-strain  curves  are  all  presented  in  the  figure. 
It  is  seen  that  the  use  of  a  flat  or  rising  stress-strain  curve  does  not  have  any  significant  effect 
on  the  midpoint  displacement  responses.  The  applied  load  has  been  non-dimensionilized  with 
respect  to  the  limit  load  PL',  which  is  given  by 


Pl 


4  T  Bb2csy ' 

jpr-J 


(4.1) 


where  B  is  the  thickness,  b  is  the  untracked  ligament  (b=W-a),  oy  is  the  yield  strength  and  S  is  the 
span  between  the  supports.  The  displacement  contour  at  the  last  load  step  for  the  rising  curve 
mode  is  presented  in  Figure  A1  of  Appendix  A. 


Plastic  Zone  Development 

The  shape  of  the  plastic  zone  was  similar  for  both  the  flat  and  rising  stress-strain  models. 
However,  the  stress  and  strain  levels  beyond  the  25%  strain  level  were  quite  different  for  the  two 
cases  because  of  the  different  load  paths.  Figure  4.2  shows  the  development  of  the  von  Mises 
stresses  and  equivalent  plastic  strains  in  the  DT-350WT  material  with  the  rising  curve  model,  at 
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different  load  levels.  For  the  sake  of  completeness,  the  corresponding  results  from  the  flat  stress- 
strain  curve  model  are  presented  in  Figure  A2  of  Appendix  A.  For  the  von  Mises  stress  plots,  all 
regions  with  stresses  above  the  yield  stress  of  450  MPa  are  shown  in  red,  and  the  area  that  has  not 
yielded  is  shown  in  blue.  For  the  equivalent  plastic  strain  plots,  all  other  areas  are  shown  in  blue. 
Also,  the  same  view  is  used  for  all  the  figures  so  a  direct  comparison  of  the  shape  and  size  of  the 
plastic  zone  can  be  made.  It  is  seen  that  the  shapes  and  size,  of  the  plastic  zones  using  the  von 
Mises  stress  or  equivalent  plastic  strain  are  very  similar.  Also,  as  expected,  the  plastic  zone  size 
increased  with  load.  In  addition,  the  plastic  zone  tends  to  rise  upwards  from  the  tip  of  the  crack, 
which  is  similar  to  the  experimentally  observed  behaviour  [2]. 

Figure  4.3  shows  the  plastic  zone  that  has  strained  over  15%  (when  necking  starts)  of  two 
load  levels.  In  comparison  with  the  plastic  zones  at  the  corresponding  load  levels  in  Figure  4.3, 
it  is  seen  that  necking  occurs  in  only  a  small  area  around  the  crack  tip. 

Plastic  Zone  Radius  (ry)  Stretch  Zone  Width  (SZW)  and  Ratio  r-y/SRW 

The  plastic  zone  radius,  ry,  and  the  stretch  zone  width  (SZW)  at  various  load  levels  were 
measured  from  the  plastic  zone  configurations.  The  plastic  zone  was  given  by  the  distance  from 
the  crack  tip  to  the  farthest  point  of  the  plastic  zone.  The  stretch  zone  width  (SZW)  was  given  by 
the  distance  from  the  crack  tip  to  the  interface  of  the  plastic  zone  with  the  side  of  the  crack  (see 
Figure  1.3).  Figures  4.4  and  4.5,  respectively  show  the  variation  of  the  plastic  radius  and  SZW 
with  load,  for  the  flat  and  rising  stress-strain  curve  models  of  the  DT-350WT  specimen.  The 
figures  show  that  the  plastic  zone  radii  and  stretch  zone  widths  obtained  from  the  two  models  are 
almost  identical  as  for  the  shapes  of  the  plastic  zones,  as  discussed  above.  The  ratios  of  the  plastic 
zone  radius  to  stretch  zone  width  (ry/SZW)  were  also  computed  and  presented  in  Figure  4.6. 
Again,  the  responses  from  the  two  constitutive  models  are  very  similar,  and  the  ratios  vary  from 
50  to  90.  Due  to  the  fact  that  the  configurations  of  the  plastic  zone  predicted  by  the  flat  and  rising 
curve  constitutive  models  were  similar,  only  the  rising  curve  (strain  hardening)  model  was  used 
for  the  remaining  specimens,  as  this  provided  more  realistic  strain  levels. 
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4.2  DT  Specimen  with  304  Stainless  Steel  Materials 

Midpoint-Displacement 

The  predicted  load-displacement  response  of  the  midpoint  of  the  beam  is  shown  in 
Figure  4.7.  The  displacements  of  the  DT-304SS  specimen  are  smaller  than  those  of  the  DT- 
350 WT  specimen.  For  instance,  at  1.1  PL,  the  midpoint  displacement  with  respect  to  the  end  point, 
DL,  is  approximately  equal  to  0.3mm  compared  to  0.63mm  for  the  DT-350WT  specimen  (see 
Figure  4.1),  since  the  applied  load  for  the  DT-304SS  specimen  (7.49  kN)  was  less  than  the  load 
applied  to  the  DT-350WT  specimen  (16  kN).  The  displacement  contours  at  the  last  load  increment 
are  presented  in  Figure  A3  of  Appendix  A. 

Plastic  Zone  Development 

Figure  4.8  shows  the  development  of  the  von  Mises  stresses  and  equivalent  plastic  strains 
in  the  DT-304SS  specimen.  In  the  von  Mises  stress  plots,  all  regions  with  stress  above  the  yield 
stress  of  240  MPa  are  shown  in  red  and  the  area  that  has  not  yielded  is  shown  in  blue.  In  the 
equivalent  plastic  strain  plots,  all  areas  with  plastic  strains  greater  than  0.01%  are  shown  in  red 
and  all  other  areas  are  shown  in  blue.  As  before,  the  same  view  is  used  for  the  plots  to  facilitate 
direct  comparison.  The  plastic  zone  configurations  at  various  load  levels  follow  the  same  pattern 
as  described  for  the  DT-350WT  specimens.  The  area  around  the  crack  tip  that  was  strained  by 
more  than  15%  was  insignificant.  So  in  Figure  4.9  only  the  area  around  the  cracktip  that  is  strained 
by  more  than  5%  is  shown.  It  is  seen  that  this  area  is  very  small  compared  to  the  overall  plastic 
zone  size. 

Plastic  Zone  Radius  (ry),  Stretch  Zone  Width  (SZW)  and  Ratio  ry/SZW 

The  plastic  zone  radius  and  SZW  at  various  loads  were  measured  as  described  in  Section 
4.1.  Figure  4.10  shows  the  variation  of  plastic  radius  with  load.  For  comparison,  the  plastic  zone 
radii  of  the  DT-350WT  specimen  are  also  shown.  As  shown,  the  plastic  radii  for  the  DT-304SS 
specimen  follow  the  same  trend  as  for  the  DT-350WT  specimen,  and  the  plastic  radii  for  the  two 
specimen  are  very  close,  with  those  of  the  DT-304SS  specimen  being  slightly  larger.  Figure  4.1 1 
show  the  SZW  variation  with  load.  Again,  the  trends  of  the  results  for  the  DT-304SS  and  DT- 
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350WT  specimens  are  similar.  However,  the  values  of  the  SZW  for  the  DT-304SS  specimen  are 
generally  smaller  than  those  of  the  DT-350WT  specimen.  The  variation  of  the  ratio  of  the  plastic 
radius  to  the  stretch  zone  width  is  shown  in  Figure  4.12,  where  a  plot  of  the  ratio  for  the  DT- 
350WT  specimen  is  also  shown.  Both  ratios  increase  slightly  with  load,  however,  the  range  of  the 
r/SZW  for  the  DT-304SS  specimen  is  140-1 50  compared  to  50-90  for  the  DT-350WT  specimen, 
this  is  expected,  since  the  304SS  material  was  more  ductile  than  the  350WT  steel. 

4.3  NS  Specimen  with  304  Stainless  Steel  Material 

Midpoint  Displacement 

The  predicted  load-displacement  response  of  the  midpoint  of  the  beam  is  presented  in 
Figure  4.13.  As  before,  DL  refers  to  the  displacement  of  the  midpoint  with  respect  to  the  end 
point  and  DS  the  displacement  with  respect  to  the  support.  Due  to  the  fact  that  the  load  was 
applied  as  a  point  load  oat  the  midpoint,  significant  local  deformation  of  the  midpoint  was 
observed.  To  reduce  the  local  deformation  effects  it  would  have  been  necessary  to  distribute  the 
load  over  a  few  elements  near  the  middle  of  the  beam.  Because  of  the  computational  effort 
required  to  do  another  run,  an  extrapolation  method  was  devised  to  reduce  the  local  effects  as 
illustrated  in  Figure  4.14.  With  this  approach,  the  computed  midpoint  displacement  of  the  last 
load  step  (1.044PL)  was  0.18  mm  compared  to  0.23  mm  obtained  experimentally  [4].  The 
displacement  contours  at  last  load  increment  are  presented  in  Figure  A4  of  Appendix  A. 


Plastic  Zone  Development 

The  predicted  von  Mises  stresses  and  equivalent  plastic  strains  at  various  load  levels  for 
the  NS-304SS  specimen  are  shown  in  Figure  4.15.  The  plastic  zones  are  defined  by  the  regions 
with  stresses  over  240  MPa  or  over  0.01%  plastic  strain.  It  is  observed  that  the  plastic  zones 
predicted  by  using  the  stress  or  strain  criterion  are  very  similar.  The  same  view  is  used  for  all  the 
plots  to  facilitate  direct  comparison.  Figure  4.16  shows  the  region  around  the  crack  with  plastic 
strains  greater  than  5%.  Again,  these  areas  are  quite  small  compared  to  the  corresponding  overall 
plastic  zone  sizes. 
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Plastic  Zone  Radius  (ry),  Stretch  Zone  Width  (SZW)  and  Ratio  r-y/SZW 
The  variations  of  plastic  zone  radius  and  stretch  zone  with  (SZW)  with  load,  for  the  NS- 
304SS  specimen  are  presented  in  Figures  4.10  and  4.11,  respectively.  For  comparison,  the 
corresponding  results  for  the  DT-350WT  and  DT-304SS  specimens  are  also  shown  in  the  figures. 
It  is  observed  that  the  plastic  radius  and  SZW  for  the  NS-304SS  specimen  are  generally  smaller 
than  those  of  the  DT-350WT  and  DT-304SS  specimens.  The  variation  of  the  ratio  of  plastic  radius 
to  stretch  zone  width  is  shown  in  Figure  4.12  along  with  the  corresponding  results  for  all  the  other 
specimens.  The  r/SZW  ratio  for  the  NS-304SS  specimen  fall  generally  between  those  of  the  DT- 
350WT  and  DT-304SS  specimens. 
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Load  Vs.  Displacement  (PL  =  14.88  kN) 


Displacement  (mm) 

♦  Displ.  DL  -  Rising  x  Displ.  DS  -  Rising _ Displ.  DL  -  Flat _ Displ.  DS  -  Flat 


Figure  4.1:  Midpoint  Displacement  Responses  of  DT-350WT  Specimen,  with  Flat  and  Rising 

Stress-Strain  Curve  Models 
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Figure  4.2:  von  Mises  Stress  and  Equivalent  Plastic  Strain  Development,  in  DT-350WT 
Specimen  with  Rising  Stress-Strain  Curve  Model 


(a)  von  Mises  Stress  and  (b)  Plastic  Strain  at  3.19  kN  (0.21  PL) 
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Figure  4.2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Strain  Development,  in  DT- 

350WT  Specimen  with  Rising  Stress-Strain  Curve  Model  • 


(c)  von  Mises  Stress  and  (d)  Plastic  Strain  at  7.87  kN  (0.53  PL) 
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Figure  4.2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Strain  Development,  in  DT- 
350WT  Specimen  with  Rising  Stress-Strain  Curve  Model 


(e)  von  Mises  Stress  and  (f)  Plastic  Strain  at  10.81  kN  (0.73  PL) 
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Figure  4.2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Strain  Development,  in  DT- 
350WT  Specimen  with  Rising  Stress-Strain  Curve  Model 


(g)  von  Mises  Stress  and  (h)  Plastic  Strain  at  14,24  kN  (0,96  Pi 
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Figure  4.2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Strain  Development,  in  DT- 
350WT  Specimen  with  Rising  Stress-Strain  Curve  Model 


(i)  von  Mises  Stress  and  (j)  Plastic  Strain  at  16.00  kN  (1.08  PJ 
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Figure  4.3:  Plastic  Zone  Configuration  of  DT-350WT  Specimen  Showing  Area  with  Greater 

than  15%  Plastic  Strain  at  7.87  kN  (0.53  PL) 
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Figure  4.4:  Plastic  Radius  of  DT-350WT  Specimen 
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Stretch  Zone  Width 


Load  (x  PL) 


A— 350WT  -  Rising  (PL  =  14.88kN)  — K-350WT  -  Flat  (PL  =  14.88kN) 


Figure  4.5:  Stretch  Zone  Width  of  DT-350WT  Specimen 


Ratio  of  Plastic  Radius  to  Stretch  Zone  Width 
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- 350WT  -  Flat  (PL  =  14.88kN) 


Figure  4.6:  Ratio  of  Plastic  Radius  to  Stretch  Zone  Width  of  DT-350WT  Specimen 
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Figure  4.7:  Midpoint  Displacement  Response  of  DT-304SS  Specimen,  Rising  Curve  Model 
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Figure  4.8:  von  Mises  Stress  Development  in  DT-304SS  Specimen  with  Rising  Stress-Strain 

Curve  Model 


(a)  at  3.16  kN  (0.46  PL)  and  (b)  at  5.92  kN  (0.87  PL) 
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Figure  4.8  Continued:  von  Mises  Stress  Development  in  DT-304SS  Specimen  with  Rising 

Stress-Strain  Curve  Model 
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Figure  4.8  Continued:  von  Mises  Stress  Development  in  DT-304SS  Specimen  with  Rising 

Stress-Strain  Curve  Model 


(e)  at  8.10  kN  (1.19  PL) 
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Figure  4.9:  Plastic  Zone  Configuration  of  DT-304SS  Specimen  Showing  Area  With  Greater 

Than  5%  Plastic  Strain  at  6.88  kN  (1.01  PL) 
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Figure  4.10:  Comparison  of  Plastic  Zone  Radius 
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Figure  4. 1 1 :  Comparison  of  Stretch  Zone  Width 
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Figure  4.12:  Comparison  of  Plastic  Radius  to  Stretch  Zone  Width 
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Figure  4.13:  Midpoint  Displacement  Response  of  NS-304SS  Specimen  With  Rising  Stress- 

Strain  Curve 
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Figure  4.14:  Correcting  for  Local  Deformations  in  NS-304SS  Specimen 
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Figure  4.15:  Plastic  Zone  Configurations  of  NS-304SS  Specimens  With  Rising  Stress-Strain 

Curve  Model 

(a)  von  Mises  Stress  and  (b)  Plastic  Strain  at  1.69  kN  (0.26  PL) 
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Figure  4.15  Continued:  Plastic  Zone  Configurations  ofNS-304SS  Specimens  With  Rising 

Stress-Strain  Curve  Model 

(c)  von  Mises  Stress  and  (d)  Plastic  Strain  at  JjJ  M  Jjj 
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Figure  4. 15  Continued:  Plastic  Zone  Configurations  ofNS-304SS  Specimens  With  Rising 

Stress-Strain  Curve  Model 

(e)  von  Mises  Stress  and  (f)  Plastic  Strain  at  4.85  kN  (0.74  PJ 
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Figure  4.15  Continued:  Plastic  Zone  Configurations  of  NS-304SS  Specimens  With  Rising 

Stress-Strain  Curve  Model 


(g)  von  Mises  Stress  and  (h)  Plastic  Strain  at  6.04  kN  (0.92  P^) 
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Figure  4.15  Continued:  Plastic  Zone  Configurations  of  NS-304SS  Specimens  With  Rising 

Stress-Strain  Curve  Model 

(k)  Plastic  Strain  at  7.79  kN  (1.18  PL) 
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Figure  4. 16:  Plastic  Zone  Configuration  of  NS-304SS  Specimen  Showing  Area  With  Greater 

Than  5%  Plastic  Strain  at  4.85  kN  (0.74  PJ 
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5.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Summary  and  Conclusions 

This  report  focussed  on  a  study  to  investigate  plastic  zone  development  in  three  point 
bend  fracture  specimen.  The  finite  element  methodology  was  utilized  to  provide  an 
understanding  of  the  plastic  zone  size  (radius)  to  crack  tip  blunting  (stretch  zone  width),  in  order 
to  be  able  to  predict  the  upper  limit  where  elastic  plastic  fracture  becomes  invalid.  Finite 
element  analysis  of  three  test  specimen  configurations  were  performed.  The  test  specimens  were 
(i)  a  standard  dynamic  tear  (DT)  test  specimen  with  350WT  steel;  (ii)  another  standard  DT 
specimen  with  304  stainless  steel;  and  (iii)  a  non-standard  three  point  bend  specimen  with  304 
stainless  steel.  The  ABAQUS  finite  element  code  was  used  to  perform  incremental  elastic-plastic 
analysis  of  the  specimens.  Intuition  and  experimental  observation  were  used  to  develop 
approximate  stress-strain  curves  for  the  materials  beyond  the  points  where  uniaxial  stress-strain 
data  were  available  for  the  350WT  steel  and  304  stainless  steel  materials. 

The  results  computed  for  each  analysis  included  the  midpoint  displacement; 
development  of  the  plastic  zone  around  the  crack  tip  with  progression  of  load;  the  plastic  zone 
radius  ry;  stretch  zone  width  (SZW);  and  the  ratio  of  plastic  zone  radius  to  SZW.  The  results 
were  presented  in  the  form  of  contour  plots  and  charts  for  comparison  with  experiments  being 
performed  concurrently.  It  was  found  that  the  plastic  zone  radius  ry  and  stretch  zone  width  SZW 
increased  parabolically  with  load,  up  to  1.1  PL,  where  PL  is  the  limit  load.  On  the  other  hand, 
the  variation  of  the  ratio  ry/SZW  with  load  fitted  a  linear  distribution  for  all  cases.  This  ratio 
was  found  to  vary  from  50-90  for  the  DT  specimen  with  350WT  steel;  140-150  for  the  DT 
specimen  with  304  stainless  steel  and  90  -140  for  the  non-standard  specimen  with  304  stainless 
steel. 


5.2  Recommendations 

The  experimental  results  were  not  fully  available  to  this  study,  as  they  were  being 
performed  concurrently.  It  is  recommended  that  the  results  of  this  study  be  compared  to  the 
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experimental  results  when  they  are  available.  Based  on  these  comparisons  it  might  be  necessary 
to  refine  the  stress-strain  relations  beyond  the  points  where  uniaxial  test  data  available,  to  obtain 
better  accuracy.  Once  a  good  level  of  accuracy  has  been  achieved  with  the  finite  element  model, 
it  is  further  recommended  that  the  analysis  be  applied  to  several  other  specimen  configurations 
to  enable  a  parameteric  study  to  be  performed.  Such  a  study  would  lead  to  the  generation  of 
generalized  charts  that  can  be  used  for  fracture  measurements  and  design. 
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Figure  Al:  Displacement  Contours  of  DT-350WT  Specimen  with  Rising  Stress-Strain  Curve 

Model  at  1.1  PL 


Figure  A2:  von  Mises  Stress  and  Equivalent  Plastic  Zone  Development  in  DT-350WT 
Specimen  with  Flat  Stress-Strain  Curve  Model 

(a)  von  Mises  Stress  and  (b)  Plastic  Strain  at  3.19  kN  (0.21  PL) 
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Figure  A2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Zone  Development  in  DT- 
350WT  Specimen  with  Flat  Stress-Strain  Curve  Model 


(c)  von  Mises  Stress  and  (d)  Plastic  Strain  at  7.87  kN  (0.53  PJ 
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Figure  A2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Zone  Development  in  DT- 

350WT  Specimen  with  Flat  Stress-Strain  Curve  Model  • 


(e)  von  Mises  Stress  and  (f)  Plastic  Strain  at  10.81  kN  (0.73  PL) 
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Figure  A2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Zone  Development  in  DT- 
350WT  Specimen  with  Flat  Stress-Strain  Curve  Model 


(g)  von  Mises  Stress  and  (h)  Plastic  Strain  at  14.24  kN  (0.96  PL) 


98-033_rpt.doc 


ttepl  mc4347.  1*1  24*«000 
Hite*  Sliet* 


>  4  50e *02 

<  4.50e*02 

<  J.75e*G2 

<  3.00e*Q2 

<  2.258*02 

<  1 . 50a  *  02 

<  ?.52e* 01 

<  2.S7a"Ql 


max  *  8  85a*0£ 
min  *  2  8?a-01 


(i) 


stepl  »nc4947.  l»1.2<«*008 
Plastio  Sttdift  Equivalent 


>  1 ,000-04 

<  1  008-04 

<  3.338-05 

<  8.678-05 

<  S.Q0«-05 

<  3.33«-0S 

<  1  67® -05 
c  0.090*00 


max  «  1  6  6  e  *  0 1 
min  *  0  008*00 


(j) 


Figure  A2  Continued:  von  Mises  Stress  and  Equivalent  Plastic  Zone  Development  in  DT- 
350WT  Specimen  with  Flat  Stress-Strain  Curve  Model 

(i)  von  Mises  Stress  and  (j)  Plastic  Strain  at  16.00  kN  (1.08  PL) 
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Figure  A3:  Displacement  Contours  of  DT-304SS  Specimen  with  Rising  Stress  Strain  Curve 

Model  at  1.044  PL 


Figure  A4:  Displacement  Contours  of  NS-304SS  Specimen  with  Rising  Stress-Strain  Curve 

Model  at  1.044  PL 
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